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ABSTRACT 

The H I Parkes All Sky Survey (HIPASS) galaxy catalogue is cross-correlated with known low 
redshift, low column density (iVni < 10 15 cm -2 ) Lyman-a absorbers from the literature. The 
redshift-space correlation is found to be similar in strength to HIPASS galaxy self-clustering 
(correlation length SQ Mg = 6 ± 4 and so.gg = 3.1 ± 0.5 Mpc respectively). In real- 
space the cross-correlation is stronger than the galaxy auto-correlation (correlation length 
ro.ag = 7.2 ±1.4 and ro. gg = 3.5 ± 0.7 /i^q Mpc respectively) on scales from 1 — 10 h^ 
Mpc, ruling out the mini-halo model for the confinement Lya absorbers at the 99 percent 
confidence level. Provided that the cause of the strong cross-correlation is purely gravitational, 
the ratio of correlation lengths suggest that absorbers are embedded in dark matter haloes 
with masses log(M/M Q ) = 14.2 h^ , similar to those of galaxy groups. The flattening of the 
cross-correlation at separations less than ~ 600 h^Q kpc could correspond to the thickness of 
filaments in which absorbers are embedded. This work provides indirect statistical evidence 
for the notion that galaxy groups and large-scale filaments, particularly those that comprise 
gas-rich galaxies, are the dominant environments of low column density Lya absorbers at 
z = 0. 

Key words: intergalactic medium, quasars: absorption lines, galaxies: statistics, large-scale 
structure of universe 



1 INTRODUCTION 

At the present epoch less than 10 per ce nt of baryoni c mat - 
ter is located in galaxies (Cole et al. 2001; Zwaan et al. 2003). 
The majority of baryons therefore lie in the intergalactic medium 
(IGM), and knowledge of the distribution of these baryons is es- 
sential to our understanding of how matter is organised in the 
Universe. Lyman-a (Lya) absorption towards distant quasars re- 
veals the presence of neutral hydrogen (Hi) in the IGM. At low 
redshift, observations show that these low column density (TYhi 
< 10 15 cm -2 ) absorbers exist in a varie ty of environments, in- 
cludin g large-scale filam ents fe.g. lLe Brun et alJl996t|Penton et all 
20021 [Rosenberg et alJl2003l) . galaxy groups (e.g. |L^nzetta_et_alJ 
1996t iBowen et alJ l2002|). and even voids iMcLinetalJ 2002) 
and underdense regions TGrogin & Gelleill 19981) . A* positive con- 
nection between abso rbers and galaxies has been e stablished by 
many studies l Morris et al. 1993; Lanzetta et al. 1995; Bowen et al. 
| l996tlTripD et alll998tlChen et alll998ll200ltllmpev et alll999t 
lOrtiz-Gil et alJ ll999). however most conclude that a one-to-one 
correspo ndence between a bsor bers and ind ividual galaxies does 
not exist. lCote et alJ 120051) and lPutman et alJ i200 6) provide com- 
pelling evidence that absorbers trace the cosmic web, finding that 
absorbers at large radii do not follow galactic rotation curves. Thus 
the absorbers appear to lie in the same overdensities as the galaxies, 
but not in individual galaxy haloes. 



email: eryan@ast.cam.ac.uk 



Despite the absorbers' positive association with highly biased, 
strongly clustered objects, i.e. galaxies and large-scale structures, 
1ow-A*hi Lya absorbers are among the weakest self-clustering 
objects in the low redshift Universe. The two-point correlation 
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namic simulations also measure weak clustering of small gas over- 
densities in velocity space l Dave et al. 2003), these overdensities, 
0.5 < log(pfr / Ph) < 1-5, correspond to Lya absorbers with 
13 < log(A>Hi/cm" 2 ) < 14 JPave et alJll999n . When compared 
with galaxies identified in the simulations, a filamentary structure 
emerges, where absorbers of increasing column density arise in 
closer proximity to galaxies. 

The cross-correlation function of absorbers and galaxies can 
be used to measure the extent to which the two populations of 
objects are associated. At intermediate redshifts the strong cross- 
correlation between Lyman break galaxies and C IV absorbers is 
used to argue that these two populations are in fact the same ob- 
ject lAdelberger et alj EoOSI The spatial co-incidence of C IV ab- 
sorbers and Lyman break galaxies however can equally be inter- 
preted as pre-enrichment from dwarf galaxies born from the col- 
lapse of ~ 2<7 fluctuatio ns at 6 < z < 12 (Mada u et ai]|200lt 
IPorciani & Madaul Eo05). The same degeneracy has the potential 
to plague the interpretation of the cross-correlation of absorbers 
and galaxies at z — 0. How can we tell the difference between a 
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strong cross-correlation signal due to absorbers arising in galaxy 
haloes and that due to a general overdensity of matter that sur- 
rounds galaxies? This issue is related to the claim that 10 14 < 
Nm < 10 18 cm -2 absorbers arise within the haloes of galaxies, 
based on the anti-correlation between impact parameter and Lya 
absor ber equivalent width lLanzetta et al]ll995tlchen et alJll998t 
|200jl). Solid observational evidence '|Bo ^n^rall200?lc6^^t^l1 
1200a). numerical sim ulations foave et ai^?999T^ and analytic mod- 
els (Linet_al. 2000) have now resolved this degeneracy, showing 
that the anti-correlation between absorber-galaxy separation and 
density of absorbing material is a natural consequence of the over- 
density of matter that surrounds galaxies, extending well beyond 
galaxy haloes. Since we know that absorbers and galaxies are 
largely not the same object, their cross-correlation signal can be 
used to measure the extent to which they are associated. The cross- 
correlation signal can also, subject to interpretation, provide an in- 
direct method to measure the bias of the absorbers, i.e. the rela- 
tionship between Hi column density and underlying dark matter 
density. 

In the Press-Schechter (PS. IPress & Schechtedl 19741) formal- 
ism, the ratio of the bias of two populations of object s is equal to 
the ra t io of t he cross- to auto-correlation function (e.g. Mo & White 
1 19961 12002I) . 1 Therefore, by (i) knowing the characteristic total 
(dark plus baryonic) halo mass of a galaxy population, (ii) mea- 
suring the galaxy auto-correlation function, and (iii) measuring the 
absorber-galaxy cross-correlation function, the mass of the haloes 
in which the absorbers are embedded can be i nferred . This method 
has been successfully used bv iBouche et alJ 120041) . They cross- 
correlated Mg II absorbers with Luminous Red Galaxies (LRGs) to 
show that the Mg II absorbers are embedded in haloes with masses 
~ 2 — 8 x 10 11 M Q , consistent with the expectation that ~ 40 
per cent of Mg II absorbers arise in Damped Lya (DLA) systems 
]Rao et ai]l2006h . and that DLAs are expected to have total halo 
masses of that order. DLAs are gravitationally collapsed objects 
and at low redshift are consist ent with the local g alaxy population 
weighted by Hi cross section iZwaan et al. 2005). Lower column 
density absorbers are not collapsed and may not follow the same 
bias trends. At intermediate redshifts the distribution of Lya ab- 
sorbe rs reveal structures on scales up to 17 h^ \ Mpc iLis ke et alJ 
2000), at these redshifts the bi as and clus tering properti es of 
Lya a bsorbers are well modelled ICen et alJl998tlMcDonald et alJ 
2002). However the physics of the IGM increases in complexity 
towards z = 0. As noted by iDave et ail 120031) the correlation 
function of gas overdensity no longer follows the correlation func- 
tion of H I optical depth at z = 0, showing that the bias is not 
as predictable as at higher redshifts. Some models however place 
absorbers in individual mini-haloes, which do have a well defined 
bias. 

Mini-halo models predict that absorbers arise in grav itation- 
ally confined ga s within dark m atter mini-haloes (Ikeuchi 1986; 
lReeslll986t Imo & Morrislll99l) . Since mini-haloes are less mas- 
sive than galaxy haloes, the absorber-galaxy cross-correlation is 
expected to be weaker than the galaxy auto-correlation. The mini- 
halo model is not an ideal description of Lya absorbers. Too many 
haloes per unit redshift are required to account for the observed 
density of absorption lines, since the mini-haloes have small spa- 
tial cross sections. 



1 In this paper the function that compares the same population of objects 
is referred to as the auto-correlation, and the function that compares two 
populations is referred to as the cross-correlation. 



The absorber-galaxy cross-correlation function is instrumen- 
tal in establishing whether absorbers are associated with individ- 
ual mini-haloes, galaxies or large-scale structure. Although other 
statistics have been used by other authors, such as nearest neigh- 
bour and galaxy density distributions, the properties of the cross- 
correlation function can be related to the underlying dark matter 
density. The cross-correlation is also a clean statistic as it does not 
rely on assumptions about correlated pairs. Mo rris et alJ al 993 ) is 
the only z = absorber-galaxy cross-correlation function pub- 
lished to-date that includes a full three-dimensional calculation. 
They used 17 absorbers along the 3C 273 line-of-sight and found 
that the absorber-galaxy cross-correlation is weaker than the galaxy 
auto-c orrelation on scales from 1 — 10 h$Q Mpc. iMo & M orris 
( 1994) explain this finding with a combination of 20 — 30 per cent 
of absorbers arising in galaxies and the remainder arising in mini- 
haloes. 

In this paper the cross-correlation function is calculated us- 
ing 129 absorbers along 27 lines-of-sight from the literature and 
5,3 17 galaxies from a blind 21-cm emission-line survey. The galaxy 
and absorber data are described in Section [2] The details of the 
cross-correlation and halo-mass calculation methods are described 
in Section|3| The subsequent results are outlined in Section|4] in- 
cluding the redshift- and real-space correlation functions. The na- 
ture of the Lya absorber correlation signal and its relationship to 
galaxies, groups of galaxies and large-scale filaments is discussed 
in Section [5] with a final summary given in Section |6| To com- 
pare with correlation functions in the literature Ho = 100 km s~ 
Mpc -1 is used throughout. 



2 DATA AND RANDOM SAMPLES 
2.1 HIPASS Galaxies 

The H I Parkes All Sky Survey (HIPASS) is a blind survey for ex- 
tragalactic 21-cm emission, thus it provides a census of gas-rich 
galaxies at 2 = 0. Extragalactic sources in HIPASS have been cat- 
alogued in the southern (5< +2°, M eyer et al.,[2004) and northern 
(2°<<5<+25.5°, Wong et al. 120051) parts of the survey. The com- 
bined HIPASS catalogues cover 29,579 square degrees and the he- 
liocentric velocity range from 300 to 12,700 km s _1 , and contain 
a total of 5,317 galaxies. HIPASS has a velocity resolution of 18 
kms -1 and opti cal ma tching shows a lcr-positional accuracy of 
1.4' (Doyle et al. l2005t) . corresponding to less than 4 h~[ Q \ kpc at 
the minimum absorber distance. The HIPASS galaxy catalogue is 
available at <www . aus-vo . org>. 

The uniform coverage of HIPASS over such a large fraction 
of the sky makes it an attractive survey for a cross-correlation 
function, as no corrections for gaps in the sur vey are re quired and 
edge effects are minimised. A recent study bv lCooper et al] 1 2005) 
finds that edges and holes in a survey can have a strong effect on 
measurements of galaxy environment. HIPASS also identifies gas- 
rich dwarf and low surface brightness galaxies, often missed by 
magnitude-limited optical surveys. In addition, the dark matter halo 
masses of HIPASS galaxies have been modelled as a function of 
H I-mass JMo etall2 005). which makes for a straightforward anal- 
ysis of the bias and clustering properties of HIPASS galaxies. 

Random galaxy samples were produced from a Poisson sky 
distribution with random velocities weighted by the HIPASS se- 
lection function (see Figure 0. The random samples each contain 
2 x 10 4 sources. 
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Figure 1. Top panel: Velocity histogram of galaxies in the combined 
HIPASS catalogues, with the HIPASS selection function overdrawn. Bot- 
tom panel: Histogram of velocities from the Lyct absorber literature sample. 



selection function should be constructed from the sum of redshift 
path lengths of the 27 lines-of-sight, but this information is not 
available for all lines-of-sight. Most (91/12 9) of the absorbers in th e 
sample are taken from lPenton et alJ fcOOfJ) or lPenton et alj 12004) . 
The path length available for the detection o f Lya absorbers for 
a superset of these sight-lines (see figure 6 of Pentoi Tet alj|2004|) 
is almost a factor of two higher at low redshift (z = 0.01) com- 
pared with that at the HIPASS redshift limit (z = 0.04). This trend 
is reflected in the shape of the histogram of Lya absorber veloc- 
ities in Figure Q Since the full absorber sample used in this pa- 
per is a compilation of data from different observing programs and 
instruments as well as data reduction, line selection and measure- 
ment techniques, weighting the velocity distribution of the random 
sample will minimise systematic effects of the non-uniform veloc- 
ity distribution. The final random absorber samples contain around 
5,000 absorbers each. Results using the different random absorber 
samples show no significant variation within the quoted uncertain- 
ties. 



2.2 Lya absorbers 

Low redshift, low column density Lya absorbers were se- 
lecte d from the S p ace Telescope Im aging Spectrograph (STIS, 
iBowen etai]l2002l: IPenton et ail l2004lwilliger et al l l2006t) and 
the Goddard High Resol ution Spectrograph (GHRS. Ilmnev et alJ 
ll999tlPenton et alj|20 00) observations in the literature. These ab- 
sorbers were chosen to overlap with the combined HIPASS cata- 
logues, i.e. <5<+25° and heliocentric velocity between 1,000 and 
12,700 km s _1 . The lower limit of 1,000 km s" 1 was applied as 
many lines-of-sight have significant Galactic absorption features. 
Only l ines with s ignificance level great er than 4.5ct were in cluded 
from llmpev et alJ il999h and 4<r from IPenton et alJ l2000l) . con- 
sistent with what these authors considered to be solid detections. 
The column densities calculated with Doppler parameter b = 30 
km s _1 are adopted from Penton et al. 1 2000 , 2004). The column 
densities of the absorbers in the literature sample lie in the range 
12.4 1 ^ loRfAWcm " 2 ) ^ 14.81, although the 10 absorbers 
from llmpev et alJ 1 19991) only have equivalent width measurements 
(in the range 0.176 to 1.809 A). For consistency, each line-of-sight 
in the sample contains absorbers from a single reference only (the 
most recent reference is used). There are 129 absorbers from 27 
lines-of-sight in the full STIS plus GHRS literature sample. 

Although a number of low redshift Lya absorbers have been 
detected with the lower resolution Faint Object Spectrograph (FOS, 
|jannuzietaTlll998t iBechtold et"ai]|2002l) , their equivalent width 
distribution and velocity uncertainty is quite different from ab- 
sorbers detected with STIS and GHRS. Furthermore, only 11 ad- 
ditional absorption lines overlap with the HIPASS survey region, 
thus only absorbers detected with STIS and GHRS have been in- 
cluded in the sample. 

Random absorber samples were produced by generating ran- 
dom lines-of-sight on the sky populated with absorbers, where the 
number distribution of absorbers per random line-of-sight matched 
the mean (4.8) and standard deviation (3.0) of the absorber sam- 
ple from the literature. Lines-of-sight with zero absorbers were not 
included in either the data or random samples. The random ab- 
sorber velocities are weighted by the velocity distribution of the 
data absorbers. No evolution is expected or observed in the ab- 
sorber population over the HIPASS redshift interval 0<z<0.04 
IPenton et afll2004 . The data absorber sample however does not 
have a flat velocity distribution (see Figure^. Ideally the absorber 



3 CROSS CORRELATION METHOD 

The cross-correl ation function, £(<r, tt), is calculated from the 
Davis & Peebles 1 1983) estimator, 



AG(a, tt) n R c 

4(0", 7T) = TT - ^ - / r 1, 



(1) 



RG(a, 7t) uag 

where AG (a, tt) is the number of data absorber-galaxy pairs with 
projected separation, a, between a — 5a/ 2 and a + 8a/ 2, and ra- 
dial separation, tt, between, tt — Sn/2 and tt + Sn/2. RG(a,n) 
is the number of pairs consisting of a random absorber and a data 
galaxy. The projected separation between two objects, Oij, is con- 
verted from angular separation, 9ij, by Oij = (veL+velj)/7Jo x 
tan(#ij /2). The separation along the line-of-sight is simply the ve- 
locity difference, -nij — vel;— velj \/Hq. Both 8a and Sir are set 
at 0.1 Ii~[q Mpc. The sampling of 0.1 h^ Mpc is used to ensure a 
data point exists at reasonably small separations in the correlation 
function. Equation Q is calculated in the range to 100 hi 00 Mpc 
in the a and tt directions. The function is normalised by the number 
of random pairs, Urg, and data pairs, hag- Only pairs with sep- 
arations less than 90° are included in the calculation (corresponds 
to 20 7i]~ q Mpc for the lowest velocity absorbers). 2 The galaxy 
auto-correlation function is calculated in a similar way, using the 
random galaxy sample. The spherical average of £(<r, tt) gives the 
redshift-space correlation function, £(s), where s = V a 2 + tt 2 . 

The separations used to calculate the cross-correlation func- 
tion do not take peculiar velocities into account. To reconcile the 
increased velocity dispersion along the line-of-sight (the finger-of- 
god effect) the cross-correlation function is integrated along the ve- 
locity axis to produce a projected correlation function, 3(cr). 



H(<x) 



£(<7, 7r)d7T 



(2) 



In practice the upper limit of the integral is 50 Mpc, and 
other upper limits are investigated. The projected and real-space 
(i.e. distortion free) correlation functions are related via 

E{a) _ 2 f x rt(r)dr 



(3) 



2 Increasing the maximum separation to 180° produces indistinguishable 
results for the projected auto- cross-correlations for a < 10 h^ Q0 Mpc, 
variations at larger a are well within the quoted uncertainties. 
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lDavis&Peeblesll983l) . If a power law functional form is assumed 
for the real-space correlation function, £ = (r/ro) _7r , then the 
parameters ro and "/ r can be derived analytically from the power- 
law coefficients of the projected correlation function, E(er) /a — 
A(r ,^ r )o~ lr , where 

j Davis & Peeblesll983l) . At large (linear) separations . iMo & White] 
( 1996) have used the PS formalism to show that the real-space auto- 
correlation function of galaxies is related to the dark matter corre- 
lation function, £dm, via the bias, b, 

£ 99 = b\M g )i DM . (5) 

This relation has been tested using N-body simulations. An equiv- 
alent relation holds for the cross-corre lation function of two popu- 
lations of objects (e.g. IMo et al J f79931. in this case absorbers and 
galaxies, 

? a9 = b(M a )b(M B )Z DM . (6) 

Thus, the amplitude of cross- to auto- correlation function can be 
expressed as the ratio of the bias of absorber dark matter haloes, 
b(M a ), to galaxy dark matter haloes, b(M g ). 

U = (r^_Y' 9a = b SMA m 
i ag W 9 J KM g ) 

As outlined in M o & White! l2002t) . the ellipsoidal collapse model 
of bias from lSheth et alj l200ll) is more accurate than the original 
spherical model, hence it shall be adopted here. To directly compare 
the amplitudes of the cross- and auto-correlation functions, their 
slopes, "/r, ag and 7r, 99 , must be equal. To ensure this the coeffi- 
cients, ro, 99 and 7r, 99 , are calculated first, then ro, ag is determined 
by fixing 7r,a 9 =7r, 99 and using a Levenberg-Marquardt non-linear 
least squares fit. For the auto-correlation function the uncertainties 
are calculated using jackknife resampling, where the variance, cr|, 
is measured by calculating the correlation function N (24) times, 
with all data and random points from 1 hour of RA removed each 
time. 

i 

Uncertainties in the cross- and absorber auto- correlation functions 
are also calculated using jackknife resampling, with one line-of- 
sight removed each time. Jackknife resampling takes into account 
cosmic variance within the HIPASS volume and represents the un- 
certainties more accurately than yN errors, which are not indepen- 
dent. 

In summary, by simply measuring the ratio of the cross- to 
auto-correlation functions, together with known halo biases and 
halo masses of galaxies, equation Q can be used to calculate the 
mass of haloes in which Lya absorbers are embedded. This method 
relies on the assumption that the relative bias can infer a mass for 
the absorber halo. 



4 RESULTS 

The auto- and cross- correlation functions, £(er, 7r) 99 and 
£(<r, vr) a9 , are given in Figures [2] and [3] The finger-of-god effect 
due to increased velocity dispersion along the line-of-sight is ap- 
parent for the galaxy auto-correlation, but is much more exagger- 
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projected separation, a (h l00 ~ 1 Mpc) 
Figure 2. The galaxy auto-correlation function. The function is calculated 
with a resolution of 5a = 8n = 0. 1 h^[ 00 Mpc and is smoothed with a box- 
car width of 9 pixels in this diagram. The smoothing level of 9 pixels is used 
to emphasize features on > 1 h^Q Mpc scales. Moderate changes to the 
sampling and smoothing levels do not change the these features. Jackknife 
realisations show typical pixel variations of the order £ 99 < 0.2. 

£(CT,7T)ag 




projected separation, a (h 100 ' Mpc) 
Figure 3. The absorber-galaxy cross-correlation function. Smoothing as per 
Figure |2| Jackknife realisations show typical pixel variations of the order 

iag < 0.5. 
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Figure 4. The absorber-galaxy cross-correlation function (§(s)aa, circles 
& solid line) and the galaxy auto-correlation function (§(s)go, crosses & 
dashed line) in redshift-space, s. The five numbers indicate the number of 
absorber-galaxy pairs that contribute to the first five £(s) a g points in the 
diagram. See text for power-law fit description. Note that the y-axis range 
varies from figure to figure. 



ated for the absorber-galaxy cross-correlation. 3 The same effect is 
seen in absorber-galaxy cross-correlation functions from simula- 
tions JPave et alJll999l) . and is interpreted as the draining of gas 
from low-density regions into collapsed structures. Figures[2|and[3] 
also reveal that at very small separations the auto-correlation signal 
is stronger, however, the cross-correlation maintains a higher value 
over a larger section of a — n space. This qualitative comparison 
suggests that absorbers are not as highly clustered as galaxies on 
< 1 /i^oo Mpc scales, though not avoiding the highest density re- 
gions (which are occupied by galaxies), the normalised chance of 
finding an absorber at larger scales (> 2 h^ 00 Mpc) is more likely 
than finding another galaxy. 

The redshift-space correlation function, £(s), is calculated 
from the spherical average of £(<r, tt). Both the auto- and cross- 
correlation functions are plotted in Figure |4] Power-law fits 
to the functions £(s) fls = (s/so,oo)~ 7s ' 99 and £(s) a g = 
{s/so,ag) ' s ' a3 are made over the range 1 — 10 h^Q Mpc for 
consistency with the projected correlation function fits (see be- 



low). The function parameters are sn, ss 

J S ,gg = 1.4±0.5, SO, OS 



±4 fojQQ Mpc, a nd 7. 



3.1 ± 0.5 



h 100 M P C > 



0. 8±0.2 . 

These results can be compared directly with Morris et alj ll993l) . 
Their pure Hubble flow galaxy auto-correlation (figure 7a) is 
marginally higher in amplitude compared with the HIPASS galaxy 
auto-correlation function. This is expected due to the weaker clus- 
tering strength of H I-selected galaxies and is discussed in much de- 
tail in lMeverN2003l) . 4 The absorber-galaxy cross-correlation func- 

3 The elongation along the velocity axis is not dominated by a particular 
line-of- sight, nor the fact that absorbers are detected along lines-of-sight, 
nor the absorbers' n(z). These possible effects have been investigated by 
calculating the cross-correlation along each line-of-sight separately, calcu- 
lating the cross-correlation with one absorber from each line-of-sight only, 
and considering only absorbers and galaxies with velocities less than 5,000 
kms . In each of these test cases the elongation along the velocity axis 
remains. 

4 Mey er uses only the southern HIPASS catalogue, the lLandv & Szalavl 
1 1993) estimator and an integration upper limit of 25 h~ nn Mpc, whereas 
both southern and northern HIPASS catalogues, the Davis & Peebles Hl983) 
estimator and an integration upper limit of 50 hTn Mpc are use here. The 
Meyer results lie within the uncertainties quoted here. 



Figure 5. The projected absorber-galaxy cross-correlation function (circles 
& solid line) and the galaxy auto-correlation function (crosses & dashed 
line). The five numbers indicate the number of absorber-galaxy pairs that 
contribute to the first five H(cr) a9 points in the diagram. See text for power- 
law fit description. 
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Figure 6. Projected cross- and auto- correlation functions compared at dif- 
ferent u values as a function of integration limit, tt, see equation J3- 



tion of lMorris et all [l 993 ) based on 17 absorbers along the 3C273 
line-of-sight is much weaker than the redshift-space correlation 
presented here, based on 129 absorbers along 27 lines-of-sight. The 
large uncertainty in the cross-correlation amplitude (so, a g = 6 ± 4) 
reflects the large variation from one line-of-sight to another, and 
emphasises the need for multiple lines-of-sight in this type of anal- 
ysis. 

As described in the method, integrating the correlation func- 
tion along the velocity axis removes the effect of peculiar veloci- 
ties along the line-of-sight and allows parameters of the real-space 
correlation function to be calculated. It is clear from Figure [3] 
that collapsing £(<j, tt) along the tt axis will produce a stronger 
signal in the absorber-galaxy cross-correlation at projected sepa- 
rations greater than 1 h^Q Mpc. The resulting auto- and cross- 
correlations are given in Figure [5] and show that the absorber- 
galaxy cross-correlation is indeed stronger than the HIPASS galaxy 
auto-correlation for the projected spatial range 1 — 10 hi 00 Mpc. 

To test the robustness of the projected correlation function, the 
upper integration limit in Equation ^2j, tt, is varied between and 
50 

^100 Mpc. Both the auto- and cross- correlation functions are 
relatively insensitive to the choice of integration limit. Lowering 
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Figure 7. The line-of-sight absorber-galaxy cross-correlation function (cir- 
cles) and the galaxy auto-correlation function (crosses). 



the upper limit of the integral does lower the amplitude of both cor- 
relation functions slightly, however H(cr) a9 remains greater than 
H(cr) 99 over a large area of parameter space. This trend is apparent 
in Figure [6] where the solid line (cross-correlation) remains above 
the dashed line (auto-correlation) for displayed values of a (dif- 
ferent shades). Galaxy-galaxy pairs with radial separations greater 
than 50 h^Q Mpc account for only 6 per cent of the total number 
of data pairs in the auto-correlation calculation. 

Parameters of the projected and real-space correlation func- 
tions are calculated by the method described in section [3] The 
power-law fit is only made over the range 1 — 10 h^Q Mpc since 
the cross-correlation flattens significantly at small separations and 
the uncertainties are high at large separations. A correlation length 
of fo.gg = 3.5 ± 0.7 /ij~ Q Mpc and slope of j r = 1.9 ± 0.3 is 
found for the auto-correlation, and ro, a9 = 7.2 ± 1.4 hi 00 Mpc 
for the cross-correlation. These parameters are plotted on the figure 
as H(cr)/(j = A(r ,~/ r )a for both the auto- and cross- cor- 
relation functions, fixing 7r,a S = 7r for reasons described in the 
method. Using equation these parameters yield a bias ratio of 
4.0 ± 1.3. The geometric mean Hi mass of galaxies contributing 
to pairs in the range 1 ^ a ^ 10 h^g Mpc is log(A/Hi/M©) = 8.8 
hioo' w ' tn a geometric standard deviation of 1.1. This corresponds 
to a halo mass of \og(M/MQ) = 1 1.0 /ij~ 00 with 1-cr lower and up- 
per li mits of log(M/M Q ) = 10.8 and 11.6 respectively jMo et alj 
2005). Combining this mean galaxy halo mass, toget her with the 
bias ratio calculated above and the linear bias relation I Mo & White 



2002) gives a value of log(M/M ) = 14.2 (with 1-cr lower 
and upper limits of 13.6 and 14.5 respectively) for the mass of 
haloes in which Lyo? absorbers are embedded. It should be noted 
that log(MHi/Mo) = 8.8 h^ is lower than the geometric mean 
Mhi of all HIPASS galaxies because the cross-correlation func- 
tion is dominated by galaxies located nearby, which on average are 
lower in H I mass. Nevertheless the bias relation at z — is a rela- 
tively flat function of halo mass for halo masses less than 10 12 Mq, 
so the choice of M g has little effect on the final result. 

The integral of £(<r, 7r) along the a axis is included for com- 
parison with single line-of-sight absorber-galaxy cross-correlation 
functions that are expressed in terms of velocity separation. An in- 
tegration limit of a — 50 h^ Mpc is used. This function is given 
in Figure the large error bars reflect the significant variations 
from one line-of-sight to another in the already weakened correla- 
tion signal alon g the veloc i ty axi s due to the finger-of-god effect. 
As noted in IWilliger et alJ feOOd) this calculation produces auto- 



Figure 8. The line-of-sight absorber-galaxy cross-correlation function (cir- 
cles) and the galaxy auto-correlation function (crosses) integrated to cr = 1 
hi 00 Mpc only to compare with the cylindrical-style cross-correlation func- 
tions from Chen et al. (2005) and Williger et al. (2005). 
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Figure 9. Absorber auto-correlation function (circles) with Penton et al. 
(2004) two-point correlation function, £(Aij), for comparison. 



and cross- correlation functions of similar amplitudes. Direct com- 
pa risons can be mad e wit h the cylindrical-style £(Au) functions 
of IChen et alj <2005l) and lWillieer et alJ 120061) by integrating to 
cr = l h^QQ Mpc (although Williger et al. have a projected separa- 
tion limit that increases from 1 to 1.6 h^ Mpc along the redshift 
cylinder) producing the function, 



B(,r) = 



a — IMpc 



£(<x', n)da', 



(9) 



given in Figure |S| Interestingly the cross-correlation drops below 
the auto-correlation fu nctio n for a < 1 ftTnn M pc in agreement 
with lChen et"ail fc005l) and lWillieer et alj 12006), however within 
the error bars, most points overlap. 

Although the mean number of absorbers per line-of-sight is 
only 4.8, it is still possible to measure the absorber auto-correlation 
in individual lines-of-sight, 



£00 = 



AA(ty) n AR 
AR(-k) riAA 



(10) 



where AA{tt) is the number of data absorber pairs with separation, 
■k, and AR(tv) is the number of random absorber pairs. Random 
absorbers are generated as described in Section l2~2l Lines-of-sight 
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Figure 10. The projected absorber-galaxy cross-correlation function for ab- 
sorbers with log (AWcm -2 ) > 13.24 and log (AWcm) ^ 13.24. The 
cross-correlation from the full absorber data set and galaxy auto-correlation 
functions from Figure 151 are also plotted for comparison (circles with solid 
line and crosses with dashed line respectively). 

with a single absorber are not included in the calculation. Figure|9| 
shows that abso rbers correlate wit h themselves ver y weakly, in 
agreem ent withlPenton et alj l2004h . shown, and also llmpev et alJ 
J 19991) and lWilliser et alJ (2006). This agreement is anticipated as 
these absorbers are drawn from those data sets, but still forms a 
good consistency check as these absorbers represent the first few in 
each line-of-sight, so there is less opportunity to form pairs. 

To test the relationship between absorber column density and 
correlation strength with galaxies, the absorber sample is split into 
two subsets defined by the me dian colum n density, log (i Vm/cm - ) 
= 13.24. For consistency with lPenton et alJ<2002tl2004 a Doppler 
parameter of 30 kms" 1 is adopted for absorption lines from 
llmpev et alJll999l) to convert the quoted equivalent width measure- 
ments into column densities. Random absorber catalogues were re- 
generated for these sub-samples based on the new velocity distri- 
bution, and new mean and standard deviation of each line-of-sight. 
Figure [Tol gives the cross-correlation for absorbers above and be- 
low the median iVm and shows definite evidence for segregation 
of lower and higher column density systems. The cross-correlation 
for absorbers with log (iVm/cm" 2 ) > 13.24 is somewhat stronger 
than that from the full sample of absorber at projected separations 
less than 1 hy^ Mpc, although the uncertainties overlap. Low col- 
umn density absorbers (log (iVm/cm -2 ) ^ 13.24) show a weaker 
correlation with galaxies, however the uncertainties are large and 
the results are consistent with no correlation on scales less than 1 

Kob M P C - 



5 DISCUSSION 

The locations of 1ow-7Vhi Lya absorbers and gas-rich galaxies have 
been cross-correlated at z = 0, giving a redshift-space clustering 
signal similar in strength to that of galaxy self-clustering. The pro- 
jected cross-correlation removes the increased velocity dispersion 
along the line-of-sight and from it the real-space clustering can be 
derived. The cross-correlation function is a powerful statistic as it 
can be related to the underlying dark matter distribution. Here the 
absorber cross-correlation results are discussed in relation to mini- 
haloes, galaxy groups, large-scale filaments, numerical simulations 
and the environs of gas-rich galaxies. 



5.1 Ruling out mini-haloes 

A typical mini-halo has a circular velocity, v c , of 30 km s _1 (e.g. 
lMo&Morri3T9 94). which corresponds to a halo mass and bias of 
6.8 x 10 9 M© and 0.73 respectively. Using Equation Q and the 
auto-correlation result, such haloes would have a cross-correlation 
length of 3.3 h^g Mpc. Comparing this correlation length with 
the absorber cross-correlation length (ro. ag = 7.2 ± 1.4 h^ 
Mpc) rules out mini-haloes for the confinement of low-iVm Lya 
absorbers at z = at the 99 percent confidence limit. 



5.2 Comparison with galaxy groups 

The ratio of real-space clustering strengths inferred from the pro- 
jected correlation functions imply that Lya absorbers are embed- 
ded in haloes with masses in the range 13.6 < log(M/M0) < 14.5 
/iT/qo . This halo mass range is of the same order as the median dy- 
namic al mass of g alaxy groups in HIPASS, log(M/M0) ~ 13.8 
/iT/qq lStevensl2005l) . This similarity warrants a careful comparison 
between correlation properties of absorbers and groups. 

The cross-correlation between groups and galaxies has been 
meas ured with both t he Sloan and 2dF galaxy redshift sur- 
veys lYang et ail 120051) . The most massive groups (halo mass 
logCM/M©) ^ 13.8 /iT/qq) display the most elongation along the 
7r axis due to larger velocity dispersions, similar to that seen in the 
absorber cross-correlation. The shape of the group-galaxy cross- 
correlation resembles tha t of the absorber-g alaxy cross-correlation 
at a < 2 and n < 10 <Yang et alJl2005l figure 1 - all groups 
and faintest galaxies), although the group-galaxy correlation level 
is higher. No 2-D flattening signature of peculiar motions due to 
infall is seen in the absorber-galaxy cross-correlation, whereas this 
signature is apparent for groups. The projected correlation is simi- 
lar in slope at 1 < a < 10, but the flattening of the absorber-galaxy 
correlation, leads to a much stronger group-galaxy correlation at 
o < 1. 

The main difference between the correlation properties of ab- 
sorbers and groups is their self-clustering. Absorber self-clustering 
(see Figure [5) appears to be weak. The weak clustering could be 
due to the measurement technique, since absorber self-clustering 
measurements are restricted to the line-of-sight axis and the clus- 
tering signal is elongated, and therefore diluted along this axis. The 
cross-correlation measured along the redshift axis has almost three 
times the amplitude as that measured perpendicular to the line- 
of-sight (compare Figures [5] and 0. Interestingly if the absorber 
auto-correlation (Figure [5) were increased by a factor of three it 
would be comparable to the redshift-space galaxy auto-correlation 
(Figure|4j- N-body simulations and analytic models can be used to 
measure absorber self-clustering in real-space. At least in the case 
of the warm-hot intergalactic medium (WHIM, slightly higher col- 
umn densities than the absorbers considered in this paper) the real- 
space cl ustering o f absorbers at r — 1 h^ Mp c could be a factor 
two less IValageas et alT2 002 ) or two more iDave et all200ll) than 
that of galaxies. 

Groups of galaxies on the other hand are strongly self- 
clustered, but not necessarily as strong as their constituent galax- 
ies. More massive haloes are more strongly clustered and con- 
tain more galaxies per halo. Thus there are many more galaxies 
than haloes in strongly clustered environments, and a comparable 
number of galaxies and haloes in the field. Therefore, purely by 
relative numbers, galaxies in strongly clustered environments will 
dominate their correlation function, whereas the correlation signal 
of groups will be more evenly weighted. This weighting can lead 
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to the self-clustering signal of groups being weaker than that of 
their constituent galaxies. Naturally, group self-clustering depends 
on the mass of haloes considered. Padilla et al. 12004 show that 
the redshift-space correlation length of groups identified in the 2dF 
galaxy redshift survey increases in amplitude from so = 5.5 /ij~ o 
Mpc for the entire sample, to 12.6 h^ Q Mpc for groups with a me- 
dian mass of 1.1 x 10 14 /ij~ Q Mr (compare d with so = 6.82 /ij~ g 
Mpc for 2dF galaxies from Hawkins et al. 2003). The projected 
auto-correlation of this larger m ass group sample is quite simi- 
lar in amplitude and shape to the lYang et alj 12005) group cross- 
correlation function described above. Hence it appears that the self- 
clustering properties of groups are likely to be strong, but could be 
weaker than galaxies. On the other hand, the self-clustering proper- 
ties of Lya absorbers are likely to be weak, but could be as strong 
as galaxies. 

Absorbers, galaxies and groups each trace the underlying dark 
matter halo distribution in different ways. By comparing the clus- 
tering properties of absorbers to that of other objects, the dark mat- 
ter haloes in which absorber reside can be revealed. As found by 
many authors (see Section 0, the possibility that Lya absorbers 
arise in individual galaxy haloes is unlikely. At a distance of 20 
^iqn M P C > HIPASS is sensitive to galaxies with Mm > 2 x 10 s 
^ioo M®. Considering only galaxies and absorbers within this dis- 
tance, the nearest galaxy to 50 per cent of absorbers is greater than 
1 h^Q Mpc in redshift-space, clearly not within the halo of an indi- 
vidual galaxy. Galaxy groups give a £(<r, n) diagram that is similar 
in shape to the absorber-galaxy cross-correlation. In addition, the 
ratio of real-space clustering lengths imply that Lya absorbers are 
embedded in haloes with masses similar to that of galaxy groups. 
These similarities provide some evidence for the notion that Lya 
absorbers are associated with galaxy group environments. 

5.3 Nature of the correlation signal compared with 
simulations 

Taking the projected absorber-galaxy cross-correlation (Figure [5J 
at face value, it can be concluded that stepping 1 to 10 h^ Q Mpc 
away from an H I-selected galaxy, the normalised chance of find- 
ing an absorber is more likely than finding another galaxy. This 
likelihood is in addition to the fact that Lya absorbers are ~40 
times more frequent than galaxies at zero redshift. 5 The exact cov- 
ering factor of neutral gas 1 h^ ( f Mpc away from a galaxy how- 
ever requires a value for the size of the absorbers, or an assessment 
of non-detections. At smaller separations the cross-correlation flat- 
tens with respect to the galaxy auto-correlation, indicating that at 
< 1 h^ ( l Mpc the normalised chance of encountering a 1ow-A?hi 
Lya absorber is less than that of another galaxy - perhaps giving 
way to higher column density systems. Flattening of the absorber 
auto-correlation function at a few hundred kiloparsecs is also seen 
in analytic models and N-body simulations of WHIM absorbers 
JPave et aliboOltlValageas et alJl2002l) . This length (~ 300 kpc) 
corresponds to the characteristic thickness of the WHIM filaments. 
A fairer comparison can be made with £(cr, n) diagrams for the 
cr oss-c orrelation of galaxies with low column density absorbers 
in |PaveetalJil999l) . Their results suggest the correlation function 
turns over at small separations, but the interpretation is limited by 
the resolution of the simulations. The results presented here show 

5 The number dens ity of abs orbers with 13.1 < log (Nm/cm~ 2 ) < 
14.0 is dN/dz=1.85 iPenton et al. 2004), the value for gas-rich galaxies is 
dN/dz=0.045 IZwaan et all2005l) . 



that the cross-correlation flattens at separations less than 600 h^ 00 
kpc. It is possible that this separation corresponds to a characteristic 
thickness of the large-scale structures in which low column density 
Lya absorbers are embedded. 

Modelling the low redshift Universe is a complicated busi- 
ness. Ideally the results presented here should be compared with 
directly with £(<j, n) and from simulations and models at zero 
redshift. Future simulations may be able to distinguish the signa- 
ture of absorbers embedded in different large-scale structures, for 
example that of groups or filaments. The essential effect that needs 
to be tested is whether the elongation along the velocity axis seen 
in the cross-correlation is due to absorbers that are embedded in 
a large gravitational potential. Another important insight offered 
by simulations and models is measuring absorber self-clustering in 
real-space. Given the large elongation of the absorber-galaxy cross- 
correlation in redshift-space, and the limitation that absorber self- 
clustering currently can only be measured along the line-of-sight, 
it is possible that the self-clustering strength of Lya absorbers has 
been understated. 

5.4 Connection with gas-rich galaxies 

Although the HIPASS galaxies are simply used as test parti- 
cles to represent the underlying density of matter, since they are 
gas-rich galaxies, ga s-rich environments are deliberately selected. 
Stevens et alj 120041) have shown that the H I content of galaxies 
in groups detected in HIPASS is no different to field galaxies, 
in contrast w ith optic ally selecte d groups, which are Hi-deficient 
IVerdes-Montenegro et aljl200lh . Thus the intragroup medium of 
H I-selected groups is expected to have a lower ionization fraction, 
leading to a high er num ber of Lya absorbers. This may explain 
why llmpev et alj 1 19991) find that absorbers in the vicinity of the 
Virgo cluster lie preferentially in regions of intermediate (rather 
than high) galaxy density, since the outskirts of galaxy clusters 
are known to be in excess of H I-rich g alaxies by a factor of three 
(Waugh et al. l2002l) . lstocke et alJi2005T) find that only 60 per cent 
of Lya absorbers with log^m/cm -2 ) > 13.5 (corresponding to 
detected Ovi absorption lines) arise in optically selected galaxy 
groups. Although the difference between optical and H I-selected 
group s in th i s case may be due to different sensitivity limits as 
IStocke et alj J2005I) consider galaxy groups with more than one 
L* galaxy, whereas HIPASS is sensitive to much fainter galaxies 
(Doyle et al. l2005h . 

All galaxies in a representative sample of HIPASS galax- 
ies (Meurer et al. 2005) have been detected in Ha emission, and 
are therefore undergoing star formation. These results emphasize 
that interstellar hydrogen is a co-requisite for star formation and 
that HIPASS gala xies can be cate gorised as emission line galax- 
ies. Interestingly, IChen et all 1200 5) finds a stronger Lya cross- 
correlation for emission line galaxies, compared with absorption 
line galaxies, a result that suggests the environs of gas-rich galax- 
ies are preferred by Lya absorbers. 

Despite the large number of absorbers and lines-of-sight, this 
study is limited by the shallow redshift coverage of HIPASS. Jack- 
knife errors take into account cosmic variance within the HIPASS 
volume, but not outside the volume. Infrared galaxy counts show 
that the local area (~200 h^Q Mpc in linear extent, to a distance of 
~150 /Ijqq Mpc) around the south Galactic pole is underdense by 
30 per cent iFrith et all2003l) . lLin et al]l2000h suggest a minimum 
redshift path length, Az, of 10 for statistics that discriminate be- 
tween absorbers arising in galactic haloes and those arising in the 
comic web. The literature sample used in this paper has Az ~ 1 
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only. HIPASS galaxies also have different clustering properties to 
optically selected galaxies <MeveJ2003l) 

For all these reasons it would certainly be worth repeating 
this study with other galaxy surveys. Future observations designed 
to measure absorber-galaxy clustering to the same extent as seen 
here would require a survey of galaxies to projected distances of 10 
h^QQ Mpc, i.e. for nearby absorbers with velocities around 10,000 
kms -1 , a radius of almost 6 degrees would be required around 
each quasar, with a depth of at least 100 h^Q Mpc along the line- 
of-sight. 



6 SUMMARY 

In this paper 129 low column density (Nm < 10 15 cm -2 ) Lya 
absorbers along 27 lines-of-sight have been cross-correlated with 
5,317 gas-rich galaxies at z = 0. A positive association be- 
tween absorbers and galaxies is found to separations of at least 
10 hi 00 Mpc in redshift space. The redshift-space auto- and cross- 
correlation functions are found to be similar in strength, with cor- 
relation lengths so, ss = 3.1 ± 0.5 and so, afl = 6 ± 4 h^ nn Mpc re- 
spectively. These results are contrary to that of M orris et alj j 19931) 
based on 17 absorbers along a single line-of-sight, which found a 
weaker cross-correlation (the uncertainty on so, ag is large however, 
so a weaker correlation for absorbers can only be ruled out at the 53 
percent confidence limit). The cross-correlation is found to depend 
on the column density of the absorbers: higher column density sys- 
tems (log (A^Hi/cm -2 ) > 13.24) are more strongly correlated with 
galaxies, especially at projected separations less than 1 h^ Mpc. 

On a one-to-one basis, out to a distance of 20 h^ Mpc, 
50 per cent of Lya absorbers are located within 1 h^Q Mpc of 
a galaxy, and the other half are located 1 — 3 h^ 00 Mpc from 
the nearest galaxy. These results agree with many othe r stud- 
ies iBowen et alll99dl2002t|Tripp et alll998l : llmpev et alll999l 
IPenton et alJl2002t ICote et alj|20 05^ that the majority of Lya ab- 
sorbers do not arise in the haloes of individual galaxies. Bearing 
this result in mind, the cross-correlation function can then be used 
to measure the extent to which absorbers and galaxies are related. 

The amplitude of the absorber-galaxy cross-correlation de- 
pends on how it is integrated. The absorber-galaxy £(<r, tt) dia- 
gram is very elongated along the line-of-sight. Integrating along the 
projected separation (a) axis produces auto- and cross-correlations 
with comparable amplitudes. Restricting this integral to 1 h^Q Mpc 
leads to a cross-correlation that is slightly weaker than the auto- 
correlation at small separati ons. Thus results from cylindrical-style 
cross-correlation functions iChen et alfcool Iwillieer et al] 2006 ) 
must be interpreted with caution. 

A more conventional approach (as employed in galaxy clus- 
tering analysis) is to integrate along the line-of-sight (n) axis to 
remove the increased velocity dispersion and from it derive the 
real-space clustering, and underlying dark matter distribution. The 
absorber-galaxy cross-correlation has been measured in real-space 
for the first time here. A correlation length of ro, as = 7.2 ± 1.4 
^ioo is found, compared with ro, 9£ , = 3.5 ± 0.7 h^[ 00 Mpc 
and slope of 7 r = 1.9 ± 0.3 for the galaxy auto-correlation. Mini- 
halo models on the other hand predict a weaker cross-correlation. 
The interpretation of these results relies on the assumption that the 
increased velocity dispersion seen in the cross-correlation is due 
to a larger gravitational potential. These results rule out mini-halo 
models for the confinement of low-iVm Lya absorbers at z — at 
the 99 percent confidence limit. The ratio of real-space clustering 
strengths, together with the bias relation (Equation imply that 



Lya absorbers are embedded in haloes with masses in the range 
13.6 < \o%{MIMq) < 14.5 h^Q , similar to that of galaxy groups. 
The cross-correlation of massive groups with galaxies produces 
a £(<7, 7r) diagram similar in shape to that of the absorber-galaxy 
cross-correlation. Thus on average Lya absorbers tend to follow 
the same cosmic web of material that embeds both galaxies and 
groups. This picture is supp orted by careful analysis of individual 
systems iBowen et all2002l) . It is also possible that absorbers arise 
in filaments of large-scale structure, although the cross-correlation 
amplitude of filaments has not been measured using observations 
nor predicted from models or simulations. 

The flattening of the cross-correlation function at projected 
separation less than ~ 600 h^ ( j kpc could correspond to a charac- 
teristic thickness of the filaments in which Lya absorber are embed- 
ded. A large discrepancy exists between individual Lya absorber 
sizes derived from photo-ionization models using metal lines asso- 
ciated with each absorber (10— 70 pc. lRiebv et alj2002llTriDP et alJ 
2002) and their minimum transverse sizes derived from quasar line- 
of-sight pair analysis (of order 300 /i7/ g kpc. IPenton et alJ l2002t 
iRosenberg et all2 0031. In agreement with Rosenben^UiiTl2"oo"3l) . 
this discrepancy could be reconciled if the minimum transverse size 
corresponds to the thickness of th e filament rat her than the indi- 
vidual absorber size. Furthermore, IPenton et alj 120021) identify a 
filament consisting of eight absorbers and five galaxies at least 20 
h^l Mpc across and 1 Mpc thick. This example fits with the ~ 
600 /i^oo kpc filament thickness inferred from the cross-correlation 
function. 

The cross-correlation results presented in this paper provide 
indirect evidence that Lya absorbers arise preferentially in gas-rich 
galaxy groups and filaments. Over half the galaxies in HIPASS oc- 
cur in groups lStevensl2 005) and gas w ith loK(iVHi/cm~ 2 ) > 13.2 
covers 50 per cent of galaxy filaments l Stocke et al. 2005) - a pic- 
ture consistent with the ubiquitous nature of Lya absorbers. 
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